\Vcd\ from D Meson Leptonic Decays 

Heechang Na, 1 Christine T. H. Davies, 2 Eduardo Follana, 3 G. Peter Lepage, 4 and Junko Shigemitsu 5 

(HPQCD Collaboration) 

1 Argonne Leadership Computing Facility, ANL, Argonne, IL 60439, USA. 
2 SUPA, School of Physics & Astronomy, University of Glasgow, Glasgow, G12 8QQ, UK 

3 Departamento de Fisica Teorica, Universidad de Zaragoza, E-50009 Zaragoza, Spain 
4 Laboratory of Elementary Particle Physics, Cornell University, Ithaca, NY 14853, USA 
Department of Physics, The Ohio State University, Columbus, OH 43210, USA 

(Dated: March 8, 2013) 

We present an update of the D meson decay constant fo using the Highly Improved Staggered 
Quark (HISQ) action for valence charm and light quarks on MILC Nf =2 + 1 lattices. The new 
determination incorporates HPQCD's improved scale r 1 * 2+1 = 0.3133(23)fm, accurately retuned 
bare charm quark masses and data from an ensemble that is more chiral than in our previous 
calculations. We find fu = 208.3(3.4)MeV. Combining the new fo with D —¥ fiv^ branching 
fraction data from CLEO-c, we extract the CKM matrix element |V ct j| = 0.223(10) e:E p. (4)j a t. • This 
value is in excellent agreement with \V c d\ from D semileptonic decays and from neutrino scattering 
experiments and has comparable total errors. We determine the ratio between semileptonic form 
factor and decay constant and find [/^ >_>,r (0) //u], t = 3.20(15)GeV _1 to be compared with the 

experimental value of [/+ "^(O) / /o] = 3.19(18)GeV _1 . Finally, we mention recent preliminary 
but already more accurate D — > [iv^ branching fraction measurements from BES III and discuss 
their impact on precision \V c d\ determinations in the future. 

PACS numbers: 12.38. Gc, 13.20.He 



I. INTRODUCTION 

Determinations of individual elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix allows for many 
cross checks and consistency tests of the Standard Model. 
In most cases there are several processes that can be used 
to extract the same CKM matrix clement each involving 
very different experimental and theory inputs. For the 
CKM matrix element \V cd \, PDG2010 [lj quotes values 
coming from D — > tt,Iv semileptonic decays and from 
neutrino/antineutrino scattering. The HPQCD collabo- 
ration recently published a new calculation of \V c d\ that 
reduced errors in the semileptonic decay determination 
by more than a factor of two Q, making it competitive 
with the neutrino scattering result. In the current arti- 
cle we present a third, independent determination based 
this time on D meson leptonic decays. We find a value 
for \V c d\ in complete agreement with the other two deter- 
minations and with comparable total errors. 

The branching fraction for the leptonic decay of a 
charged D or D s meson via a virtual W boson is given 
to lowest order by, 
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where m; is the charged lepton mass and q — d or s. 
Electromagnetic corrections to this formula are known 
and routinely taken into account by experimentalists in 
their analyses [3, |j| . Equation (QJ tells us that determi- 
nation of | V c d | from D leptonic decays requires theory to 
provide the D meson decay constant fjj which is a pure 



QCD nonperturbative quantity. The first lattice QCD 
calculations of fn and fo, that included sea quark contri- 
butions were carried out by the Fermilab Lattice & MILC 
collaborations @ and this predated experimental studies 
of these decays. Subsequent experimental measurements 
were consistent with the lattice predictions within errors 
that were more substantial then than they are today for 
both theory and experiment. The initial lattice calcula- 
tions employed an effective theory approach (the heavy 
clover action Q) for the charm quark on the lattice. In 
2007 the HPQCD collaboration introduced the Highly 
Improved Staggered Quark (HISQ) action which repre- 
sents not only an extremely accurate lattice quark action 
for light quark physics, but also serves as an accurate 
relativistic action for heavier quarks Q. The HISQ ac- 
tion has since been used very successfully in simulations 
involving the charm quark such as for charmonium, and 
for D and D s meson decay constants and semileptonic 
form factors [a&G3- In reference HPQCD published 
the first f n , fx, Jd and Jd 3 results from HISQ valence 
quarks, including HISQ charm quarks, on the MILC Asq- 

all with sub 2% errors 



Tad Nf — 2 + 1 lattices [11 



At 



around the same time experimental measurements of D 
and D s meson leptonic decay branching fractions were 
improving significantly |L2h14| |. And by the middle of 
2008 we were facing an interesting situation where there 
was good agreement between experiment and theory for 
fo but a close to 4c discrepancy in fr> s . Further im- 
provements and scrutiny became crucial. 

The largest systematic error for fr> s in reference Q 
came from the uncertainty in the scale r\. HPQCD 
was using an r\ extracted from T splittings namely 
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TABLE I: Simulation details on three "coarse" and three 
"fine" MILC ensembles. 



Set 


n/a 


m,i/m s (sea) 


N con f 


Ntsrc 


I? x N t 


CI 


2.647 


0.005/0.050 


1200 


2 


24 3 x 64 


C2 


2.618 


0.010/0.050 


1200 


2 


20 3 x 64 


C3 


2.644 


0.020/0.050 


600 


2 


20 3 x 64 


FO 


3.695 


0.0031/0.031 


600 


4 


40 3 x 96 


Fl 


3.699 


0.0062/0.031 


1200 


4 


28 3 x 96 


F2 


3.712 


0.0124/0.031 


600 


4 


28 3 x 96 



ri = 0.321(5) fm with 1.56% errors [15]. In 2010 
HPQCD published a much more accurate r± determina- 
tion, r\ — 0.3133(23), based on several physical quanti- 
ties and an improved continuum extrapolation (from 5 
lattice spacings) [l6| . A change in the scale affects quan- 
tities such as fr> s in two ways: 1. the bare strange and 
charm quark masses must be retuned on each ensemble 
and 2. the conversion from dimensionless decay constant 
(e.g. in units of r{) to the decay constant in physical 
units is modified. In reference Q HPQCD updated its 
value for fjj a together with and fx using the new 
r\. Although and fx hardly shifted at all upon go- 
ing from old to new r\ , the updated fo 3 came out about 
2.3cr (3%) higher than before. As a consequence, taking 
into account also that experimental results were changing 
and moving closer to theory numbers, the discrepancy in 
fu s between theory and experiment has now shrunk to 
a 1.6(7 effect. Reference [9[ did not present a new calcu- 
lation of fo ■ Instead we took the previous ratio Jd s / Id 
from reference @ and combined this with the new fjj a 
to estimate a new fu- 

In this article we complete the process of switching to 
the new r\ scale for meson decay constants and present 
a direct calculation of fjj consistently using the new 
scale. Since the time of reference Q experimental er- 
rors in the D — > /i, branching fraction have improved 
from - 7.8% down to ~ 4.3% in the case of CLEO-c [ij] 
and new even more accurate measurements are appear- 
ing now from BES III [l?} • Together with the new fjj of 
this article with its ~ 1.66% error, one can now extract 
a \Vcd\ from D meson leptonic decays that is as accurate 
as those from semileptonic decays or neutrino scattering 
and that promises to become even more precise in the 
near future. 



II. THE LATTICE SETUP 

Tabic I lists the three coarse (a « 0.12fm) and three 
fine (a « 0.09fm) MILC ensembles used in this study 
together with some lattice details. And in Table II we 
show the values for valence quark masses. For fjj we 
have focused more on ensuring better control over chi- 
ral extrapolations by adding a more chiral fine ensemble 
(Set F0) rather than going to finer lattices as we did 



TABLE II: Valence quark masses 



Set 


ami 


am a 


arric 


CI 


0.0070 


0.0489 


0.6207 


C2 


0.0123 


0.0492 


0.6300 


C3 


0.0246 


0.0491 


0.6235 


F0 


0.00339 


0.0339 


0.4130 


Fl 


0.00674 


0.0337 


0.4130 


F2 


0.0135 


0.0336 


0.4120 



in reference [9| for fjj s . The bare charm quark mass is 
tuned using the physical rj c mass adjusted for the absence 
of electromagnetic, charm sea and annihilation contribu- 
tions in our simulations which leads to a target value of 
M target = 2 .985(3) GeV [3 rather than the experimen- 
tal Value of M^p = 2.980(1) GeV. Most of the charm 
quark mass tuning had been done already in reference 
101 for our D — > K,lv studies. For the present calcu- 
lations we needed to add tunings only on ensemble F0. 
Figure 1. shows the tuned rj c masses for all 6 ensembles. 
The bulk of the errors shown comes from the ~ 0.1% 
uncertainty in ri/a, whereas the tiny black error bars 
represent the statistical errors on each data point. A 
similar plot for tuning of the strange quark mass via the 
r] s (fictitious) meson mass is give n in Figure 3 of refer- 
ence [l3| ■ And in references (8l4l0j| we have demonstrated 
that once quark masses have been fixed by rj c and rj s 
then masses for the D and D s mesons can be derived 
with zero adjustable parameters in good agreement with 
experiment. We do not repeat those calculations here. 
However, since we have new data for the mass difference 
AMu = Md s — Md, we summarize them in an Appendix 
and compare with AMh = Mb b — Mb in the B system 
taken from reference [191 ]. 

Having fixed the quark masses we evaluated D and D s 
correlators on each of the 6 ensembles. We use random 
wall sources with a different set for each color component 
in order to improve statistical errors. In the next section 
we describe how we extract meson decay constants from 
these correlators. 



III. CORRELATORS AND FITTING 
STRATEGIES 

The decay constant fo of a pseudoscalar meson made 
out of a charm quark and a light antiquark of mass m q 
is defined in terms of the matrix element of the heavy- 
light axial vector current — v I / g 7 M 75 v I' c between the 
hadronic vacuum and the D meson state. 

(0\A tl \D)=p lt f D (2) 

Since we employ the relativistic HISQ action for all 
quarks we are able to take advantage of PCAC, as is 
routinely done for f n and fx, and express the decay con- 
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FIG. 1: Checking the tuning of the charm quark mass to 
the Tj c meson mass. Errors on the simulation results include 
statistical (black error bars) plus errors arising from the un- 
certainty in ri/a for each ensemble. The "experimental" r\ c 
mass has been adjusted to take into account the lack of anni- 
hilation and electromagnetic effects in our lattice calculation. 



stant in terms of the pseudoscalar density PS — 4 , g 75^ / c, 

(3) 



The relevant hadronic matrix element in eq.© can be 
extracted from the D meson two-point correlator, 

c D nt (t) = ^Y,J2^D(y,mU^)\o), (4) 

x y 

where $d = PS x a 3 is the same as the pseudoscalar 
density in lattice units, and is used here also as an inter- 
polating operator for the D meson. We fit C^ nt (t) to 



the form, 

C 2 D pnt (t) = J2 |6^| 2 (e- B °* + e-^^-*)) 



N D -1 



k=0 
iV'n-1 



+ |dj°I 2 (-l)*(e- B i D * + e-Etf'Cw*-*)). 



fc=0 



(5) 



The ground state amplitude b$ is related to the matrix 
element of interest as, 



]h pp- \(0\*D\D)\ 2 \(0\PS\D)\ 2 a 3 
\ b o I = o^_„3 = ^TZ ' ( 6 ) 



2M D a 3 



2M T 



and hence, 



a.fD = 



m c + m q 

Mn 



aMr. 



(7) 



Our initial goal is to extract the amplitude = \ | as 
accurately as possible. In references [2|,llCj we found that 
fit results for two-point energies and amplitudes are im- 
proved significantly if one carries out simultaneous fits to 
two-point and three-point correlators. Three-point cor- 
relators are calculated, for instance, when one studies 
D — > TTjlv semileptonic decays. For pions at zero mo- 
mentum one has, 



C^(i,T) 



^EEE 



x y 

t r# 



(8) 



where S is the heavy-light scalar density ^f c ^q in lattice 
units. Cp™^ must be fit to the form, 



N„—l n d -i 
3 k 



(T-t) 



+ E E B^e-^e-^^i-l)^ 

3 k 
N'„-l N D -1 

+ E E C jk e-^*e-^-*\-l)* 

3 k 
N' -1 N' D -l 

+ E E ^•fee- j5 r* e - B "( T -*)(-i)*(-i)( T - i ). 

3 k 

(9) 

We will only consider the region < t < T and take 
T « N t so that any contributions from mesons prop- 
agating "around the lattice" due to periodic boundary 
conditions in time can be ignored. The same energies 
and E'® appear in ([5]) and (|9]). Doing simultane- 
ous fits to Cj^ nt and C 1 ^"* places tighter constraints on 
these energies and this helps in reducing fitting errors in 
the two-point amplitudes b® . In this way the three-point 
correlator is acting like a very complicated but effective 
smearing for the propagation of D mesons. Normally this 
would also be considered a very expensive smearing, how- 
ever we already had simulation results for on five 
out of the six ensembles in Table I from the D semilep- 
tonic project published in reference 0] so we could take 
advantage of this. It was only necessary to create new 
three-point correlator data on ensemble F0 and this only 
for zero momentum pions. 

In Fig. 2 we show some results for b$ on ensemble CI 
versus the number of exponentials from simultaneous fits 
(we set No = N' D = N n = N^) and compare with fit 
results to just C^ pnt alone. One sees the improvement 
in the fitting errors coming from the simultaneous fits. 
All our fits are done using Bayesian methods [2(j. We 
use the "sequential method" , where starting from N =2 
or 3 the output from an N - exponential fit becomes the 
initial values for the subsequent (JV+1) - exponential fit. 



4 



0.2 



0.195 



0.19 



0.185 



0.18 



O SimFit result 
O 2pt fit result 



<> 



$$$$$$ 



4 



FIG. 2: Comparison of fit results for the amplitude b$ for 
ensemble CI coming from just the two-point correlator (red 
data point to the left of figure) versus results from simultane- 
ous fits to two-point and three-point correlators. 



TABLE III: Fit Results 



Set 


aM D 


afD 


aM Ds 


a .fo 3 


fojfn 


CI 
C2 
C3 


1.1395(7) 
1.1591(7) 
1.1618(5) 


0.1370(5) 
0.1421(4) 
0.1464(3) 


1.1878(3) 
1.2014(4) 
1.1897(3) 


0.1541(3) 
0.1566(3) 
0.1552(3) 


1.1245(37) 
1.1018(32) 
1.0600(16) 


F0 
Fl 
F2 


0.8096(3) 
0.8130(3) 
0.8189(3) 


0.0943(2) 
0.0966(2) 
0.1001(2) 


0.8471(1) 
0.8471(2) 
0.8434(2) 


0.1074(1) 
0.1082(1) 
0.1076(1) 


1.1385(22) 
1.1202(21) 
1.0750(14) 



In addition to the D meson decay constant Jd we have 
also accumulated new data for fjj s by studying D s meson 
two-point correlators. Here we do not have D s semilep- 
tonic decay three-point correlator data. So, our extrac- 
tion of the relevant amplitude b® s was carried out from 
just the two-point correlators. Since statistical errors are 
smaller for D s than for D mesons, this lack of ability to 
carry out simultaneous fits in the case of D s was not a 
serious problem. In Table III we list all our fit results for 
aMu, afo, aMjy s , afn 3 and the ratio /d s //d- 



IV. CHIRAL AND CONTINUUM 
EXTRAPOLATION 

The next goal is to extrapolate the entries for fn in 
Table III to the continuum and chiral limit. The latter is 
defined as the limit m q lrn s — > 1/27.4, or using m s /m c — 
1/11.85 from reference [2lj], the limit m q /m c -> 1/(27. 4x 
11.85). We carry out the simultaneous chiral/continuum 
extrapolation using continuum partially quenched heavy 
meson chiral perturbation theory (PQHMChPT) [H-Q 



augmented by lattice spacing dependent terms. This is 
the same formalism employed recently in our Jb and fs 3 
determinations [l9j . We write, 

f D = A(l + Sf + [analytic])(l + [discret.]). (10) 

The chiral logarithm term Sf is taken from the original 
literature on PQHMChPT fp, and is also summa- 
rized in the Appendix of [19] . As in that reference we 
take, 

[analytic] = 

/3 (2m u + m s )/m c + /3im q /m c + /3 2 (m q / m c ) 2 , 

(11) 

where m u (m q ) is the sea (valence) light quark mass. rh c 
is the AsqTad charm quark mass tuned to the rj c meson 
made out of AsqTad charm quark and antiquark, and is 
the appropriate charm quark mass to use for sea quarks. 
We take rh c from reference Q where it was found that 
rh c /m c ~ 0.9 for lattices employed in the current article. 
Using ratios of bare quark masses to parameterize the 
"analytic" terms is convenient since such ratios are scale 
independent. We use the valence charm quark mass as 
the scale to measure the dominant discretization effects 
and set, 



[discret.] — co(am c ) 2 + ci(am c ) 4 



(12) 



We will call the chiral/continuum extrapolation ansatz 
given by eq.JTU]) together with (Till), © and eq.(A7) of 
reference [19( for 5f our "basic ansatz". The result of 
the extrapolation to the physical point using the basic 
ansatz is given by the green square point in Fig. 3. We 
have tested the stability of this result by modifying the 
basic ansatz in a number of ways and redoing the ex- 
trapolation. The modifications that were tried out are 
the following : 

1. dropping the term in (fTTj) 

2. adding a (m q /m c ) 3 term in (|lll) 

3. dropping the c\ term in (|12l) 

4. adding (am c ) n , n = 6, 8, 10, to IQ2} 

5. replacing a in (fl"2"|) by a x [power series in (m q /m c )] 

6. using powers of (a/fi) rather than of (am c ) in (|12[) 

7. using eq.(Al) of reference [l9| rather than (A7) for 
the chiral logarithm term Sf 

8. allowing for a 20% error in f n entering the chiral 
perturbation theory formulas 

Fig. 4 compares the extrapolation results with these mod- 
ifications in place with the basic ansatz value at the phys- 
ical point. 
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V. RESULTS 
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FIG. 3: Results for fn versus m ; /m s and at the physical point 
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FIG. 4: Testing the stability of the chiral/continuum extrap- 
olation. The numbers on the horizontal axis refer to the type 
of modification applied to the basic ansatz of (|10[1 , and 
(|12fl as listed in the text. The orange horizontal lines bracket 
the basic ansatz result, i.e. the physical point result in Fig.3. 



TABLE IV: Error budget 



Source 


fn 


fn. 


fojfn 




(%) 


(%) 


(%) 


statistics / fitting 


0.5 


0.3 


0.3 


scale n 


0.7 


0.7 




n/a 


0.04 


0.05 




continuum extrap. 


1.2 


1.2 


0.9 


chiral extrap. & Qd* Dtt 


0.7 


0.2 


0.5 


mass tunings 


0.1 


0.2 


0.2 


finite volume 


0.3 


0.1 


0.3 


Total 


1.7 


1.5 


1.1 



Table IV gives the error budget for fo 7 Jd s and 
IdJId- For all but the last two entries we use the meth- 
ods of reference [25] to isolate contributions from different 
sources that make up the total error coming out of the 
chiral/continuum extrapolations. For the finite volume 
error we take over the result from reference Q where 
an analysis was carried out comparing finite and infinite 
volume chiral perturbation theory. 

Taking all errors into account our final value for fo is, 



fn = 208.3(1.0) stai .(3.3) sl/s .MeV. 



(13) 



This is in good agreement with the previous result of 
f D = 207(4)MeV [| based on HPQCD's old n, but is 
slightly more accurate. Eq. (fT3")) represents the most pre- 
cise fn available today. 

For completeness we also give new values for fjj s and 
foj Jd, 



f Ds = 246.0(0.7W(3.5) sys .MeV, 



and 



fnjfn = 1.187(4) stat .(12; 



sys. 



(14) 



(15) 



The result for fjj s , eq. (fT4")) . is consistent with HPQCD's 
best updated value of f Ds = 248.0(2.5)MeV @ but is not 
as accurate. One sees from Table IV that the dominant 
error comes from the continuum extrapolation. In this 
respect the current calculation of /d s is not competitive 
with reference [9| which employed data from five lattice 
spacings. 

The new of eq. (fT3")) can be combined with the D — > 
/i, branching fraction from CLEO-c [l2j to extract a 
new value for \V c d\- We find, 



w< 



cd \ lepton.d 



= 0.223(10)^.(4) 



flat.- 



(16) 



The first error, which is the experimental error, domi- 
nates the total error of 4.8%. Eq. (fTfI|) agrees very well 
with HPQCD's recent determination of \V c d\ from D — > 
tt,Iv semileptonic decays 0, namely \V c d\ sem iie P .d. = 
0.225(6) e:C p.(10)i a i., where now the lattice error domi- 
nates over the one from experiment. Both leptonic and 
semileptonic determinations agree with \V c d\ = 0.230(11) 
l] coming from neutrino scattering, and all three have 
comparable total errors. 

As mentioned in the Introduction, BES III has recently 
announced preliminary results for the D —> [iis^ branch- 
ing fraction (l7j . Using their numbers we find, 

\Vcd\Zpfonl = 0.220(7W(4)iat. [preliminary]. (17) 

which agrees well with (|16p and has smaller experimental 
errors. 

Another way to check the consistency of the Standard 
Model and/or to test the lattice approach to heavy fla- 
vor physics is to consider the ratio between semileptonic 
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FIG. 5: Result for fp from this work and comparisons with 
previous work 



form factor and decay constant /+^ ir (0)//£>. We find, 
by combining eq.(l3]) with /f^O) = 0.666(29) from 
reference @, 



[f?^ w (0)/f D ] lat =3.20(15)GeV~ 1 . 



(18) 



This can be compared with the experimental ratio in 
which | V c d | cancels of [H], [26| 

[f?^(0)/f D ] exp _ =3.19(18)GeV- 1 . (19) 



Eq. (fT3]) . cq. (jT6|) and the good agreement between ([18 
and ([1"§1) arc the main results of this article. 



VI. SUMMARY 

In this article we presented a new determination of the 
CKM matrix element \V cc i\, eq. (fTB)l . made possible by an 
updated calculation of the decay constant fn, eq. (ll"3l) . 
and improved determinations of the D — » /j,, leptonic 
decay branching fraction by CLEO-c [l2j and BES III 
In Fig. 5 we compare the new fjj with HPQCD's 
previous val ue |8j a nd with results from other lattice col- 
laborations [24l |27|. l28j |. And in Fig. 6 we plot different 
results for \V c d\ including the leptonic decay determina- 
tion of this article, together with semileptonic decay and 
neutrino scattering determinations. 

In the future it will be important to continue working 
on reducing the theory errors in eq. (|18l) and the exper- 
imental errors in eQ. (|19[i . The former is dominated by 
errors in the lattice determination of ""^(O) and work 
is underway to significantly reduce them [29J. The ex- 
perimental error in eq. (|19|) comes mainly from the lep- 
tonic decay branching fraction and one can look forward 



FIG. 6: Summary of \V c j\ determinations from leptonic and 
semileptonic D decays, from neutrino scattering and from uni- 
tarity. 



to improvements there as well. In particular, the recent 
measurements by BES III ljj look very promising. The 
crucial question is whether the nice agreement seen now 
between eq. (fT8|) and eq. (fT9| will continue to hold once 
errors dip down to ^2% or below. 
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Appendix A: The D a - D Mass Difference 

In this appendix we summarize results for the mass dif- 
ference AMd = Mo s — Mjy and compare with the anal- 
ogous difference in the B system AMb = Mb . — Mb, 
where the latter was calculated in reference [19| employ- 
ing NRQCD 6-quarks. This is an interesting quantitiy 
to compare since the leading heavy quark mass depen- 
dence cancels in each of the mass differences and one is 
testing whether the subleading contributions are accu- 
rate enough to be able to distinguish between the D and 
B systems. In the difference of differences AMb — AMb 
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TABLE V: Mass Splittings in the D and B systems. The 
AMb numbers are taken from [191 ]. 



Set 


AM D [MeV] 


AM B [MeV] 


AMb - AM B [MeV] 


CI 


80.4(1.1) 


64.8(2.2) 


15.6(2.5) 


C2 


69.7(1.0) 


57.7(1.8) 


12.0(2.1) 


C3 


46.5(5) 


41.3(2.0) 


5.2(2.1) 


FO 


87.3(7) 


71.7(2.9) 


15.6(3.0) 


Fl 


79.4(7) 


61.4(2.0) 


18.0(2.1) 


F2 


57.4(4) 


47.8(1.3) 


9.6(1.4) 



are used in the B/B s and the D/D s calculations). Ta- 
ble V lists simulation results for AMb, AMb and for 
AMo — AMb- The first two quantities are plotted in 
Fig. 7 versus mi/m s . For AMb statistical errors are small 
enough so that a slight lattice spacing dependence is de- 
tected. Errors are larger for AMb and no discretization 
effects are visible. Fig. 8 shows AMb — AMb- Ones 
sees agreement with experiment at the 1 a level, with a 
a corresponding to about 5 MeV. With current levels of 
improvements to the lattice actions, 5 MeV appears to be 
the accuracy with which the HISQ action or the combined 
NRQCD/HISQ actions are able to describe charm-light 
or bottom-light boundstate dynamics. 
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FIG. 7: Comparison of the mass differences AMd = Md s 
M d and AM B = M B , - M B - 



FIG. 8: AMd — AMb, the difference of the mass differences 
in the D and B systems. 



any mistunings of the strange quark mass should also 
cancel out (identical strange and light quark propagators 
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